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ABSTRACT 


Back  flow  from  the  exhaust  of  a  chemical  laser  in  low  earth  orbit 
may  be  detrimental  to  the  integrity  and  operation  of  the  selfsame 
system.  Difficulties  arise  in  the  calculation  of  exhaust  plume  properties 
and  molecular  flux  as  the  gas  expands  from  continuum  to  free-molecular 
flow.  The  solution  of  the  governing  Boltzmann  equation  is  exceedingly 
complex;  similarly,  numerical  solutions  such  as  the  Direct  Simulation 
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I.  INTRODUCTION 


Advancing  satellite  and  spacecraft  technology  will  soon  enable  the 
deployment  of  an  array  of  space-based  laser  systems.  Among  the  lasers 
being  considered  are  chemical  lasers  of  the  hydrogen  halide  class. 
Before  their  deployment  may  be  realized,  however,  important  questions 
remain  to  be  answered.  Major  concerns  are  the  possible  contamination 
and  system  degradation  of  low  earth  orbit  spacecraft  by  laser  exhaust 
back  flow.  One  design  of  a  spacecraft  laser  is  that  of  a  long  cylinder, 
its  power  output  device  located  at  one  end,  as  in  Figure  1.1. 


Wall  ^ 
contamination 


Ambient  atmosphere 
Exhaust  gas 


Disturbance  to  laser 


Laser  beam 


King  nozzle 

Figure  1.1  Spacecraft  Laser. 

The  exhaust  gas  is  released  through  a  ring  nozzle  and  undergoes 
rapid  three-dimensional  axisymmetrical  expansion.  However,  a  number 
of  molecules  in  the  expanding  plume  may  have  thermal  velocity  compo¬ 
nents  directed  back  toward  the  spacecraft  of  sufficient  magnitude  to 
overcome  the  mean  flow  velocity.  As  certain  molecular  species  in  the 
exhaust  plume  of  a  chemical  laser  may  be  corrosive  to  the  spacecraft 
surface,  the  molecular  backscatter  from  the  jet  must  be  determined  to 
assess  potential  hazards  to  the  spacecraft. 


Considerable  difficulty  arises  in  the  analysis  of  the  flow  field  of 
the  exhaust  jet  as  it  expands  rapidly  from  continuum  flow  conditions  at 
the  nozzle  exit  plane  into  the  near-vacuum  conditions  of  low  earth 
orbit.  In  the  region  where  an  exhaust  plume  transitions  from 
continuum  to  free-molecular  flow,  a  microscopic  approach  must  be  taken 
in  the  flow  analysis.  Properly,  temporal  changes  in  molecular  distribu¬ 
tion  functions  in  physical  and  velocity  space  must  be  related  to  molec¬ 
ular  collision  parameters.  In  such  analysis  the  only  applicable  equation 
is  the  Boltzmann  equation  [Ref.  1],  viz.. 


ff,)go-  dOdu, 


(1.1) 


Here,  n  is  the  number  density,  f(u,,x,,t)  is  the  normalized  distribution 
function,  uf  is  molecular  velocity,  o  is  a  differential  collision  cross 
section,  dQ  is  an  element  of  solid  angle,  F,  is  an  external  force  per 
unit  mass,  and  du  is  an  element  of  velocity  space.  The  subscript  1 
denotes  the  molecular  collision  partner  and  the  prime  denotes  a  post¬ 
collision  value.  The  left-hand  side  of  the  Boltzmann  equation  repre¬ 
sents  the  change  in  the  number  of  particles  in  a  phase  space  element 
from  convection  of  molecules  across  the  surface  of  this  space  due  to 
molecular  velocity  u,  and  external  forces  F,  .  The  right-hand  side  of 
the  equation  represents  the  change  in  the  number  of  molecules  in  a 
phase  space  element  due  to  particle-to-particle  collisions.  The 
Boltzmann  equation  has  been  the  subject  of  extensive  investigation; 
however,  the  complexity  of  the  problem  has  thus  far  yielded  analytical 
solutions  only  to  simplified  models  of  a  specific  nature  [Refs.  2,3,4]. 

For  instance,  the  Chapman-Enskog  method  of  solution  of  the 
Boltzmann  equation  assumes  the  distribution  function  f  is  perturbed  by 
a  small  amount  from  the  Maxwellian  form.  The  distribution  function  is 
then  expressed  as  a  power  series  in  e,  a  parameter  regarded  as  a 
measure  of  mean  collision  time  [Ref.  5], 

Approximations  to  the  Boltzmann  equation  by  Hilbert,  Burnett,  and 
others  consists  of  expanding  f  in  powers  of  \/a,  where  X  is  the  mean 
free  path  at  a  point  in  the  flow  field  and  o  is  a  typical  macroscopic 
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dimension  of  the  flow  geometry.  In  the  first  approximation  to  the 
Boltzmann  equation,  the  left-hand  side  of  equation  1.1  is  omitted,  and 
the  result  is  a  nonlinear  integral  equation.  Higher  approximations 
satisfy  inhomogeneous  integral  equations  [Ref.  6].  However,  the 
validity  of  the  above  methods  is  confined  to  flow  of  small  Knudsen 
numbers,  i.e.,  o»X. 

A  numerical  approach  to  the  solution  of  the  Boltzmann  equation 
proposed  by  Bird  [Ref.  7]  is  the  Direct  Simulation  Monte  Carlo  (DSMC) 
method,  which  entails  the  computer  modeling  of  a  real  gas  flow  by 
several  thousand  simulated  molecules.  Continuum  flow  solutions  (e.g., 
the  method-of-characteristics)  are  used  to  compute  flow  field  properties 
to  the  point  at  which  continuum  flow  theory  is  no  longer  valid. 
Beyond,  flow  is  simulated  on  a  microscopic  scale,  molecule-by-molecule. 
Velocity  components  and  spatial  coordinates  of  each  of  thousands  of 
molecules  are  modified  over  small  increments  of  time  until  a  steady-state 
solution  is  determined.  A  shortcoming  of  the  Direct  Simulation  Monte 
Carlo  technique  is  that  such  a  solution  may  be  obtainable  only  after  a 
great  number  of  iterations  of  molecular  collision  algorithms,  requiring 
prodigious  amounts  of  computer  processing  time. 

In  order  to  significantly  reduce  the  computational  time  required  for 
molecular  back  flow  analysis,  while  obtaining  results  of  the  correct 
order  of  magnitude,  an  alternate  approach  to  the  problem  is  presented 
herein. 


II.  METHODOLOGY 


'A 


Consider  the  gas  jet  of  Figure  2.1,  exiting  an  axisymmetric  ring 
nozzle  into  a  near-vacuum.  The  gas  undergoes  Prandtl-Meyer  expan¬ 
sion  in  which  the  flow  is  turned,  Mach  number  increases,  and  density 
decreases  until  ambient  pressure  is  reached. 


x 


Figure  2.1  Prandtl-Meyer  Expansion. 
Region  A  of  Figure  1  is  expanded  in  this  Figure. 


At  some  point  along  any  streamline  in  such  a  flow,  there  is  a  break¬ 
down  of  continuum  flow  theory  [Ref.  8].  Although  the  transition  from 
viscous  to  free-molecular  (collisionless)  flow  is  actually  continuous,  a 
boundary  where  the  transition  begins  may  be  approximated  in  three 
dimensions  by  a  breakdown  surface.  We  wish  to  evaluate  the  magnitude 
of  molecular  efflux  from  the  breakdown  surface  toward  the  surface  of 
the  spacecraft.  In  reality,  molecules  crossing  the  breakdown  surface  in 
the  direction  of  the  spacecraft  enter  a  transition  region  in  which  molec¬ 
ular  collisions  are  highly  likely  to  occur,  as  depicted  in  Figure  2.2. 


11 


Continuum  Transition 
y|  flow  /region^ 

•Exhaust/,  7^ 


'rree  molecular 
(collisionless)  flow 


Figure  2.2  Exhaust  Flow  Breakdown. 

However,  within  the  transition  region,  translational  equilibrium  of  the 
expanding  flow  ceases  to  exist,  and  thermal  properties  of  the  gas,  in 
effect,  remain  as  they  were  at  the  onset  of  breakdown.  Therefore, 
although  molecular  collisions  occur  between  the  breakdown  surface  and 
the  free-molecular  boundary,  it  is  assumed  that  effusion  of  molecules 
toward  the  spacecraft  from  the  breakdown  surface  is  equivalent  to  that 
of  the  transition  region.  Hence,  this  analysis  is  concerned  with  deter¬ 
mination  of  flux  incident  on  the  spacecraft  from  a  defined  breakdown 
surface. 

Bird  [Ref.  9]  has  defined  a  "breakdown  parameter"  for  steady  flow 


p-u. 
pv  ds 


(2.1) 


where  v  is  collision  frequency  per  molecule,  p  is  density,  u  is  flow 
speed,  and  s  is  distance  along  a  streamline.  The  breakdown  parameter 
is  thus  a  modified  Knudsen  number  relating  mean  free  path  to  density 
scale  length.  The  onset  of  breakdown  occurs  when  the  value  of  P  is 
approximately  0.05.  We  define  r  as  the  distance  from  the  fan  vertex  in 


the  Prandtl-Meyer  expansion  and  show  that  P  may  be  expressed  explic¬ 
itly  as  a  function  of  r  and  Mach  number.  The  locus  of  the  breakdown 
surface  may  thus  be  determined.  The  nozzle  lip,  at  which  point  the 
Prandtl-Meyer  expansion  waves  are  centered,  is  a  singularity  point 
which  is  also  the  origin  of  lines  of  constant-P,  as  shown  in  Figure  2.3. 


x 


Figure  2.3  Contours  of  Constant-P. 

As  the  nozzle  lip  is  approached,  uncertainties  in  flow  field  analysis 
increase.  Here,  we  are  concerned  only  with  flow  properties  of  the  gas 
located  more  than  a  few  tens  of  mean  free  paths  from  the  nozzle  lip 
singularity.  As  the  gas  flow  proceeds  through  the  Prandtl-Meyer 
expansion,  the  value  of  r  on  the  breakdown  surface  increases  as  the 
local  Mach  number  increases.  Properties  of  the  flow  at  any  point  along 
the  breakdown  surface  may  be  determined  as  a  function  of  stream  flow 
speed.  Further,  the  molecular  flux  from  a  given  element  of  area  on  the 
breakdown  surface  is  a  function  only  of  local  Mach  number  and  the 
angle  between  the  flux  vector  and  mean  flow  direction.  Thus,  from  any 
element  on  the  breakdown  surface,  the  molecular  flux  incident  on  an 
element  of  area  on  the  spacecraft  surface  may  be  determined. 

At  some  point  along  the  breakdown  surface,  the  Mach  number  is 
sufficiently  high  that  the  backward  flux  effusion  toward  the  spacecraft 


surface  is  virtually  nil,  and  consideration  of  flux  from  points  further 
out  on  the  breakdown  surface  is  unnecessary.  This  point  is  shown  in 
Figure  2.4  at  the  intersection  of  the  breakdown  surface  defined  by 
P  =  0.05  and  the  expansion  wave  corresponding  to  this  limiting  Mach 
number. 

P  ■  .05 

yt  \ 


. .  x 

Figure  2.4  Extent  of  Significant  Portion  of  Breakdown 
Surface  as  Defined  by  the  Limiting  Mach  Number. 

The  total  incident  flux  at  any  point  on  the  spacecraft  surface  is  calcu¬ 
lated  by  integration  over  all  those  portions  of  the  breakdown  surface 
from  which  significant  flux  is  received. 


III.  SOLUTION 


A.  EXHAUST  GAS  PROPERTIES 

In  the  schematic  representation  of  a  low  earth  orbit  spacecraft 
laser  system  of  Figure  1.1,  pertinent  dimensions  shown  are  the  space¬ 
craft  radius  R1  and  half-width  of  the  exhaust  nozzle  XI.  The  exhaust 
jet  of  a  hydrogen  fluoride  laser  consists  primarily  of  four  molecular 
species:  hydrogen  fluoride,  hydrogen,  deuterium  fluoride,  and  helium. 
Laser  cavity  molar  flow  rates  and  gas  properties  [Refs.  10,11,12,13] 
assumed  for  the  calculations  which  follow  are  weighted  averages  repre¬ 
sentative  of  typical  laser  cavity  parameters  rather  than  those  of  a 
specific  operating  point.  These  parameters  are  as  follows: 

molecular  weight,  M. . 7.4  kg/kmole 

molecular  diameter,  6  .  2.52  x  10“ m 

specific  heat  ratio,  1 . 1.493 

stagnation  density,  Pj . 0.01  to  0.1  kg/m* 

stagnation  temperature,  Tj .  2200  K 

exit  Mach  number,  M„ . .  .  3.0  to  5.0 


B.  PRANDTL-MEYER  EXPANSION 


Laser  exhaust  gas  exits  the  spacecraft  from  a  centrally-located 
axisymmetric  ring  nozzle  and  undergoes  two-dimensional  isentropic 
expansion.  The  following  relations  are  assumed  valid  at  any  point  in 
the  isentropic  flow  field  [Ref.  14]: 


(3.3) 


’ 


T,  =  T[l+(r-1)MV2] 
p,  =  p[i+(r-i)Ml/2]v,(r"> 

a  =  v9RT  <35) 

Local  temperature,  density,  and  Mach  number  are  denoted  by  T,  p,  and 
M.  The  stagnation  temperature  Tt  and  specific  heat  ratio  T  are  each 
constant  in  the  flow  field.  R  is  the  gas  constant  per  unit  mass,  it  is 
the  universal  gas  constant,  and  a  is  local  sonic  speed. 

Figure  2.1  depicts  a  streamline  of  exhaust  flow,  with  Mach  number 
Ms  at  the  nozzle  exit  plane,  undergoing  Prandtl-Meyer  expansion  as  it 
leaves  the  nozzle.  The  angle  between  the  streamline  and  a  given  Mach 
wave  is  a  function  of  the  local  Mach  number,  related  by 

sin  fj,  —  1/M  (3.6) 


Also, 


p  =  e  +  n  (3.7) 

In  equation  3.7,  8  is  the  angle  between  the  spacecraft  surface  and  the 

flow  streamline;  3  is  the  angle  between  the  spacecraft  surface  and  the 

tail  of  the  Prandtl-Meyer  expansion  fan.  Owczarek  [Ref.  15]  gives  the 

total  angle  that  the  flow  is  turned  from  its  direction  at  the  exit  plane 

as 


e-e2  =  (3.8) 

where  the  angle  of  turn  is  equal  to  0,  -  02>  Values  for  v  are  given  by 
the  Prandtl-Meyer  function 

KM)  =  V(7+1)/(7H)  tan-V(7-l)(M2-1)/(7+|)  -  tan-Vj?H  (3.9) 

In  Figure  2.1,  0,=  90°  and  M,=  M«. 


_ 
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C.  THERMALLY  SCATTERED  FLUX 


it  is  assumed  that  the  molecules  in  any  small  element  of  volume  in 
the  expanding  flow  possess  a  Maxwellian  velocity  distribution.  Noeller 
[Ref.  161  shows  that  the  fraction  of  molecules  in  such  a  volume  element 
that  fly  into  the  space  angle  du,  as  depicted  in  Figure  3.1,  is 

(drtr/nT)d(Jidc  =  Tr-V!c!9-(cW-2eu  c  *)  dwdc  (3  ,0) 


Figure  3.1  Flux  Geometry. 


Here,  u  is  the  stream  velocity,  c  is  the  inertial  flow  velocity  (i.e., 
thermal  velocity  plus  stream  velocity),  and  *  is  the  angle  between  the 
direction  of  flight  and  the  mean  flow  direction.  Flow  speeds  in  this 
analysis  have  been  non-dimensionalized  with  respect  to  the  most  prob¬ 
able  molecular  thermal  speed,  cm,  where 


c„  =  n/2RT 


(3.11) 


Thus, 


u  =  Ma/V2RT  =  MVtRT/V^RT  =  v^M 


(3.12) 


•  •  \ 


►'*  **»  .*•  V'  V'V* 


a 


The  number  density  of  molecules  reaching  a  point  P  from  the  volume  t 


dnp  =  (dnT/d<y)dfl 


(3.13) 


where  dfi  is  the  solid  angle  of  the  surface  F  from  P.  Analogously,  we 
may  define  the  incident  molecular  flux  in  the  vicinity  of  P  arriving  from 
the  surface  F  as 


dip  =  (<%/dtt)dO 


(3.14) 


where 


(d|F/dcy)  =  nTcm7r-V2cse“(c2+u2“2cu  cos  ^dc  (3.15) 


Hence, 


*  /»oo 

(djp/dO)  =  I  nTcm7r-3/2c3e-(c2+u2-2cu  cos  tfdc 
J  0 

The  evaluation  of  this  integral  is  given  in  Appendix  A.  Thus, 
d|p  =  nTcm7r'3/2e”u2[(1+u2cos2^)/2+(\/7T/2)(u  cos 

(eu2cos2  ^erfc^-y  cos  ^)(u2cos2  if/  +  3/2)]  dO  (3.17) 


The  local  number  density  is  given  by 


n  =  p(  N/tt) 


where  N  is  Avogadro's  number.  From  equations  3.4  and  3.18, 


(3.18) 


p,(N/.«) 


T  [t+(r-t)MV2] 


2  /«ll/(7-l) 


(3.19) 


^  *  o  *  ’  •  *  *  v.i  *  *  '  ’  •  •  • 

•>  v  v  v  v  rv  w*  *.* 

4*  ^  *  •  a.  m  ll/J  II  Iktk  ^  A  \  k 


~  *  '  *.*  >  V  ‘ 


From  equations  3.3  and  3.11 


cm  =  V2RT,/[1+(7-1)mV2]  (3.20) 

The  total  flux  arriving  at  point  P  from  the  breakdown  surface  is  found 
by  integrating  equation  3.17  over  the  entire  solid  angle  dG  within  which 
the  surface  F  appears  from  P. 

From  equations  3.12  and  3.17-3.20,  the  flux  emanating  from  a  small 
surface  element  is  a  function  only  of  local  Mach  number  and  the  angle 
off  - the  mean  flow  direction,  </>.  The  flux  is  a  maximum  in  the  direction 
of  u,  decreasing  rapidly  as  <(>  increases.  From  equation  3.17,  values  of 
(djp/dQ)  were  calculated  for  various  Mach  numbers  as  a  function  of  1>. 
Results  are  shown  in  Figure  3.2.  Values  along  the  ordinate  are 
flux^,  M  /flux^«0.M*3'  i  •  e. ,  the  ratio  of  flux  directed  at  an  angle  ♦ 
from  the  mean  flow  direction  to  the  magnitude  of  flux  at  tp  =  0°, 
M  =  3.0.  For  example,  for  local  flow  at  Mach  3.0,  flux  at  a  back  flow 
angle  of  180  degrees  is  approximately  seven  orders  of  magnitude  less 
than  flux  in  the  stream  flow  direction.  Similarly,  at  Mach  14.5,  flux 
directly  opposing  the  stream  flow  is  nearly  eighty  orders  of  magnitude 
below  the  maximum  value. 

D.  DETERMINATION  OF  BREAKDOWN  SURFACE 

A  parameter  for  the  determination  of  continuum  flow  breakdown  is 
given  in  equation  2.1.  Falcovitz  [Ref.  17,  Appendix  B]  evaluates  the 
breakdown  parameter  P  along  streamlines  of  a  centered  rarefaction  fan. 
A  non-dimensional  gradient  is  defined  as 


whose  value  for  a  Prandtl-Meyer  fan  is 


/e(M)  =  2Vm!— 1/M(y+1) 


(3.22) 


as  a  Function  of  Mean  Flow  Mach  Number, 


Bird  [Ref.  18]  gives  the  following  relations  for  the 
speed  c,  mean  free  path  X,  and  collision  frequency  v: 

average  thermal 

c  =  2cm/v£ 

(3.23) 

X  =  (v5rr<52n)"’ 

(3.24) 

£ 

II 

(3.25) 

Substitution  of  equations  3.5,  3.11,  and  3.22-3.25  into  equation  2.1 
yields 

V(yA)(M2-1) 

P=  2(y+1)rd2n  (3.26) 

Substitution  of  equation  3.18  into  equation  3.26  gives 

V(yA)(M2-1)[l+(y-1)MV2]l/l7-,) 

2(y+l)rd2pf(N/M)  (3-27) 

To  determine  the  location  of  the  breakdown  surface,  P  is  set  equal 
to  0.05  in  equation  3.27.  We  may  now  solve  explicitly  for  the  distance 
r  of  the  breakdown  surface  from  the  vertex  of  the  expansion  fan.  As  a 
function  of  local  Mach  number, 

-  -  10V(7AXM2-1)[H(7-1)MV2l1/(T-1> 

(y+f)i2p,WM )  t3'23) 

From  evaluation  of  the  Prandtl-Meyer  function  for  the  flow  fields  under 
consideration,  it  is  found  that  the  maximum  Mach  numbers  reached  in 
the  expansion  to  the  ambient  pressure  of  low  earth  orbit  are  well  in 
excess  of  one  hundred.  However,  the  molecular  back  flow  from  areas 
along  the  breakdown  surface  with  flow  speeds  of  this  magnitude  is 
negligible.  Thus,  for  efficiency  of  computation,  it  is  necessary  to  esti¬ 
mate  the  maximum  Mach  number  Mm„ which  will  bear  significantly  in  the 
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present  analysis.  An  estimate  of  M,,,,,*  15.0  was  arrived  at  in  the 
following  manner.  For  the  gas  under  consideration,  the  Prandtl-Meyer 
function  (equation  3.9)  was  evaluated  for  M  =  15.0.  For  flow  exiting 
the  nozzle  at  Mg  =  3.0,  the  angle  through  which  the  flow  is  turned 
(equation  3.8)  is  approximately  50°.  Therefore,  the  back  flow  angle 
for  molecules  from  this  flow  to  reach  a  point  in  the  plane  parallel  to  the 
spacecraft  surface  is  approximately  40°  off  the  mean  flow  direction. 
Evaluation  of  equation  3.16  with  the  appropriate  values  gives 

log,0(flux+M0./f(ux^,0-)M „  ,s.0  ~  -30  °  <329> 

That  is,  the  flux  directed  toward  a  plane  parallel  to  the  spacecraft 
surface  from  flow  at  Mach  15.0  is  thirty  orders  of  magnitude  less  than 
the  flux  in  the  mean  flow  direction.  Moreover,  the  magnitude  of  flux 
arriving  at  the  spacecraft  surface  will  be  less  than  this  value.  Thus 
M  =  15.0  is  used  as  an  a  priori  estimate  of  M„,0II.  Subsequent  calcula¬ 
tions  confirm  that  loss  of  accuracy  is  not  incurred  using  this  estimate. 


E.  BACK  FLOW  GEOMETRY 


Molecular  back  flow  geometry  is  presented  in  Figure  3.3. 
Two-dimensional  projections  onto  the  x-y  plane  and  y-z  plane  are  shown 
in  Figures  3.4  and  3.5.  From  Phillips  [Ref.  19],  the  angle  4>  between 
the  mean  flow  stream  and  a  point  P  on  the  spacecraft  surface  is  given 
by 


COS  i/  = 


(3.30) 


where  L  is  the  vector  from  a  point  B  on  the  breakdown  surface  to  a 
point  P  on  the  spacecraft.  The  magnitude  of  u  is  equal  to  the  product 
of  local  Mach  number  and  sonic  speed.  Define  the  vector  from  the 
origin  to  P  as  P,  the  vector  from  the  origin  to  B  as  B,  and  unit 
vectors  in  the  x,  y,  and  z  directions  as  x,  y,  and  z. 
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Figure  3.5  Back  Flow  Geometry:  y-z  plane. 


From  Figures  3.3,  3.4,  and  3.5,  the  following  relations  hold: 

P  =  XX  +  Rly 

=  (XI  +  r  cos  +  [(RI  +  r  sin  £)sin  0]y 

+  [(R1  +  r  sin  0) cos  0]z 


u  -  (u  cos  0)x  +  (u  sin  8  sin  <p) y  +  (u  sin  6  cos  0) z 

a* 

L  =  P-B  =  (x  -  XI  -  r  cos  /3)x  +  [RI  -  (RI  +  r  sin  /3)sin  0]> 

~  [(RI  +  r  sin  fi)cos  0jz 


|L|=  j(x  -  XI  -  r  cos  /?)2  +  [(RI  -  (RI  +  r  sin  0)sin  0]2 

+  [(RI  +  r  sin  0)cos  0]2j 


t/2 


(3.31) 


(3.32) 


(3.33) 


(3.34) 


(3.35) 


cos  f  =  j(u  cos  0)(x  -  XI  -  r  cos  0)  +  (u  sin  5  sin  0) 

[R1  -  (R1  +  r  sin  0)sin  0]-  (u  sin  0)(R1  +  r  sin  0)cos 2  0j  /|L||u|  (3.36) 

The  solid  angle  dfl  within  which  an  element  of  the  breakdown 
surface  appears  from  the  point  P  is 


dfi  =  (  L-rW)/|L|” 


(3.37) 


where  n  is  the  unit  vector  normal  to  the  breakdown  surface  and  do  is 
an  element  of  area  on  the  breakdown  surface.  Phillips  [Ref.  20]  gives 


ndcr  =  ^  i?  drd0 

or  00 


(3.38) 


where  r  is  the  vector  along  the  Mach  wave  for  a  specified  Mach  number 
to  the  point  B  on  the  breakdown  surface  (Figure  2.4).  The  term  dr 
denotes  an  incremental  distance  along  the  breakdown  surface.  With 
reference  to  the  coordinate  system  of  Figures  3.3,  3.4,  and  3.5: 


r  =  (r  cos  0)x  +  (r  sin  0  sin  0)y  +  (r  sin  0  cos  0)2  (3. 39) 


*gjr  =  (cos  0)x  +  (sin  jS  sin  0)y  +  (sin  cos  0)2  (3.40) 


»-  ^  -v  Jfcia&aaifcg 


Thus, 

nda  =  [(-r  sin2  /S)x  +  (r  sin  /S  cos  0  sin  0)y 

+  (r  sin  j8  cos  0  cos  <£)z]  drd0 


C3.42) 


and 

L-ndcr  =  |(x  -  XI  -  r  cos  j3)(-r  sin2  /3)  +  (r  sin  0  cos  /S  sin  0) 

[R1  -  (R1  +  r  sin  /S)s?n  0]  -  (r  sin  /3  cos  /S  cos  0) 

(R1  +  r  sin  /S)cos  drd#  ^3-43) 

The  solid  angle  d$2  (equation  3.37)  may  now  be  determined,  and  the 
flux  from  an  element  of  breakdown  surface  to  an  element  of  area  near 
point  P  (equation  3.17)  calculated.  The  total  flux  is 
f'T*/2-f'rma  x 

djP  =  27T~5/2j  J  nTcme“u2  [(1+u2cos2^)/2+(Vti/2) 

rmin  rmin 

(u  cos  Tp)(eu2coa2  ^)erfc(-u  cos  ^)(u2cos2  V'  +  3/2)]  dO  (3‘44) 


Note  that  flow  symmetry  with  respect  to  the  y-axis  requires  integration 
only  from  4  =  4mln  to  4  =  90°.  Twice  this  result  gives  the  value  for 
total  flux. 


F.  SPATIAL  INTEGRATION  OF  FLUX 

The  independent  variable  used  in  calculations  thus  far  is  local 
Mach  number.  For  given  stagnation  parameters  and  exit  plane  Mach 
number,  the  continuum  flow  breakdown  surface  and  all  pertinent  flow 
properties  may  be  determined  as  has  been  shown.  The  FORTRAN 
program  'FLUX'  (Appendix  C)  was  written  to  perform  the  calculations 
presented  herein.  In  the  program,  equation  3.44  is  evaluated  over  the 
extent  of  the  breakdown  surface  contributing  flux  to  a  given  target 
area  on  the  surface  of  the  spacecraft.  Only  areas  of  the  breakdown 
surface  whose  y-components  are  greater  than  the  spacecraft  radius  R1 
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in  Figure  3.5  will  contribute  flux  to  a  fine  of  points  along  the  space¬ 
craft  surface  defined  by  y  =  R1,  2  =  0.  This  constraint  may  be 

expressed 


(R1  +  r  sin  /?)s in  <f>  >  R1  (3.45) 

Having  chosen  the  value  of  Mmo,a  maximum  value  for  distance  from 
the  nozzle  lip  rm„  is  defined  by  equation  3.28.  This  value  will  in  turn 
define  a  minimum  value  of  *  for  the  integration  procedure,  given  by 

Kin  =  sin"[RV(Rt  +  rmox)]  (3.46) 

Included  in  Appendix  C  is  a  flow  chart  of  the  program  'FLUX'. 
The  integration  algorithm  begins  with  the  input  of  constants,  gas  prop¬ 
erties,  and  spacecraft  geometry.  An  exit  Mach  number,  maximum  Mach 
number,  and  an  incremental  step  value  dM  are  chosen,  and  a  string  of 
target  points  P  are  selected  on  the  surface  of  the  spacecraft  along  the 
line  y  =  R1,  z  =  0. 

Calculations  begin  for  a  target  point  near  the  nozzle  lip.  In  order 
to  avoid  the  singularity  at  the  nozzle  lip,  the  value  of  M  is  initialized 
as  (Me*dM)  in  a  plane  near  *  =  *mlo.  The  ♦-plane  in  which  integration 
begins  is  defined  by  values  assigned  to  *m„,  and  d*.  In  this  constant-* 
plane,  Mach  number  is  incremented  by  dM  and  the  constraint  of  rela¬ 
tion  3.45  is  checked.  When  a  value  of  r  is  reached  such  that  rela¬ 
tion  3.45  is  satisfied,  the  location  of  the  point  B,  pertinent  flow 

parameters,  the  back  flow  angle  4>,  and  the  incremental  flux  from  B 
toward  the  target  point  are  calculated.  The  value  of  M  is  incremented 
and  the  process  is  iterated. 

A  flux  magnitude  of  one  molecule  per  square  meter  per  second 

arriving  in  the  vicinity  of  a  given  target  point  is  used  as  the  minimum 
value  for  which  calculations  proceed  in  a  constant-*  plane. 


That  is,  if 


(diP)M  0  <  1  per  m2  per  sec  (3.47) 

then  the  value  of  0  is  incremented  by  d#,  M  is  reset  to  M  ,  and  itera¬ 
tion  begins  in  a  new  plane  of  constant-^.  This  process  continues 
through  4  =  90°.  Total  incident  flux  near  a  given  target  point  is  equal 
to  twice  the  sum  of  the  incremental  flux  values  thus  calculated.  A  new 
target  point  is  selected  and  the  entire  procedure  is  repeated. 


IV.  RESULTS 

Calculations  were  performed  as  described  above  for  exit  Mach 
numbers  of  3.0,  4.0,  and  5.0,  and  exhaust  plume  stagnation  densities 
of  0.01  and  0.1  kg/m*.  A  constant  stagnation  temperature  of  2200  K 
was  used  in  all  calculations.  Table  1  indicates  some  of  the  limiting 
values  calculated  for  the  stated  initial  conditions.  Maximum  and 
minimum  distances  of  the  breakdown  surface  from  the  nozzle  lip  are 
given  by  rmM  and  rm,„.  Mm,„is  that  Mach  number  along  the  breakdown 
surface  at  which  the  flux  to  any  target  area  on  the  surface  becomes 
less  than  one  molecule  per  square  meter  per  second.  Note  that  the 
maximum  Mach  number  for  which  this  constraint  is  reached  is  M  =  12.8. 
Therefore,  the  a  priori  estimate  of  M  =  15.0  did  not  limit  the  accuracy 
of  back  flow  calculations.  The  two  right-hand  columns  give  the  mean 
free  path  X  at  the  nozzle  exit  plane  and  the  ratio  of  rmlla  to  X.  The 
latter  parameter  gives  a  relative  indication  of  the  "closest  point  of 
approach"  of  the  breakdown  surface  to  the  singularity  at  the  nozzle  lip. 


TABLE  1 


Exit  Plane  and 

Breakdown 

Surface 

Parameters. 

M« 

P\>  kg/m3 

rmax*  cm  rmln*crn 

^mox 

X,  cm 

rmln/^ 

3.0 

.01 

81 .53 

.1730 

11.65 

.0047 

37.1 

3.0 

.10 

13.00 

.0173 

12.80 

.0005 

37.1 

4.0 

•91 

27.63 

.5592 

9.35 

.0111 

50.2 

4.0 

.10 

3.57 

.0559 

9.85 

.0011 

50.2 

5.0 

•91 

i8.3§ 

1.4878 

8.60 

.0236 

63.0 

5.0 

.10 

2.29 

.1488 

9.00 

.0024 

63.0 

The  most  salient  data 

presented  in 

Table 

1  are  the  values 

which  in  essence  define  the  extent  of  the  breakdown  surface  contrib¬ 
uting  flux  to  the  surface  of  the  spacecraft.  Values  for  r  are  gener¬ 
ally  less  than  one  meter.  Moreover,  for  the  cases  of  higher  stagnation 
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densities  and  exit  Mach  numbers,  the  total  significant  efflux  from  the 
breakdown  surface  toward  the  spacecraft  originates  from  within  only  a 
few  centimeters  from  the  nozzle  lip. 

Figures  4.1,  4.2,  and  4.3  show  the  magnitude  of  flux  incident  on 
the  spacecraft  surface  for  exit  Mach  numbers  of  3.0,  4.0,  and  5.0, 
respectively.  In  each  figure,  values  along  the  abscissa  (distance  from 
the  nozzle  lip)  range  from  one  millimeter  to  ten  meters.  For  each  set  of 
initial  conditions  considered,  flux  in  the  plane  of  the  nozzle  exit 
reaches  a  maximum  near  the  nozzle  lip,  generally  well  within  ten  centi¬ 
meters.  The  variance  of  flux  within  this  distance  of  the  lip  is  not  more 
than  about  one  order  of  magnitude.  Beyond,  incident  flux  to  the 
spacecraft  decreases  rapidly;  at  a  point  ten  meters  from  the  nozzle  exit 
the  back  flow  is  generally  three  to  four  decades  less  than  the  maximum 
values. 

The  magnitude  of  molecular  back  flow  is  a  strong  function  of  exit 
plane  Mach  number,  exhibiting  inverse  proportionality.  The  variance  of 
back  flow  with  stagnation  density  is  less  pronounced.  Near  the  nozzle 
lip,  the  spacecraft  receives  a  slightly  greater  molecular  flux  from  the 
higher-stagnation  density  exhaust  plume.  However,  at  distances 
greater  than  one.  meter  from  the  nozzle  lip,  flux  received  by  the  space¬ 
craft  is  virtually  independent  of  stagnation  density,  and  becomes  a 
function  only  of  exit  Mach  number. 

Figure  4.4  depicts  the  sensitivity  of  flux  magnitude  to  assumed 
values  of  the  breakdown  parameter  P.  As  P  is  modified,  the  general 
locus  of  the  breakdown  surface  changes  as  shown  in  Figure  2.3. 
Maximum  flux  magnitude  appears  to  be  insensitive  to  increasing  values 
of  P.  However,  changes  in  the  shape  of  the  flux  distribution  curve  do 
occur.  As  the  assumed  value  of  P  increases,  flux  near  the  nozzle  lip  is 
increased  and  the  point  of  maximum  flux  moves  inward  toward  the 
nozzle  lip.  Beyond  the  point  of  maximum  flux,  though,  the  decrease  in 
back  flow  to  the  surface  of  the  spacecraft  decays  more  rapidly  for 
higher  values  of  P. 


Calculations  were  also  performed  to  assess  the  validity  of  approxi¬ 
mating  the  breakdown  surface  by  the  expansion  wave  corresponding  to 
M*,*,  as  shown  in  Figure  2.4.  From  results  shown  in  Figure  4.5,  it 
appears  that  little  loss  of  accuracy  is  incurred  in  making  this 
assumption. 


Figure  4.3  Molecular  Flux  in  Plane  of  Nozzle  Exit* 
Stagnation  Temperature  ■  2200  K.  Exit  Mach  Number  ■  5.0 
Plume  Stagnation  Densities  are  given  in  kg/m*. 


V.  CONCLUSIONS 


The  data  presented  in  Figures  4.1,  4.2,  and  4.3,  and  Table  1 
indicate  that  for  both  effusion  and  reception  of  molecular  flux,  the 
region  in  the  immediate  vicinity  6f  the  nozzle  lip  is  of  critical  impor¬ 
tance.  Therefore,  efforts  to  minimize  the  molecular  back  flow  to  the 
spacecraft  surface  warrant  attention  to  this  region.  In  theory,  rather 
straightforward  and  uninvolved  modifications  to  the  spacecraft  structure 
near  the  exhaust  nozzle  would  significantly  decrease  incident  back  flow. 
Comparison  of  Figures  4.1,  4.2,  and  4.3  shows  that  minimization  of 
back  flow  is  contingent  upon  raising  the  exit  Mach  number  of  the 
exhaust  jet  to  a  value  as  high  as  feasible.  With  all  other  parameters 
constant,  increasing  the  exit  Mach  number  from  3.0  to  5.0  decreases 
flux  to  any  given  area  on  the  spacecraft  surface  by  approximately  five 
orders  of  magnitude. 


Figure  5.1  Exhaust  Nozzle  Modification. 

Shown  in  Figure  5.1  is  the  addition  of  a  rim  located  at  either  edge 
of  the  nozzle,  encircling  the  spacecraft.  Benefits  of  such  an 
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arrangement  would  be  twofold.  Most  importantly,  the  effective  nozzle 
expansion  ratio  is  increased,  thereby  increasing  exit  Mach  number. 
Additionally,  the  back  flow  angle  for  molecular  flux  to  reach  a  given 
area  on  the  spacecraft  is  greater  than  that  for  the  design  in  which  the 
nozzle  exit  plane  is  flush  with  the  spacecraft  surface. 


Figure  5.2  Exhaust  Nozzle  Modification. 

Small  vertical  barriers  encircling  the  spacecraft  near  the  nozzle  lip 
as  shown  in  Figure  5.2  would  also  be  effective  in  significantly  reducing 
incident  flux.  Added  without  regard  to  increasing  area  ratio  or  exit 
Mach  number,  the  "fences"  would  act  simply  to  block  a  high  percentage 
of  the  molecular  effusion  from  the  breakdown  surface.  In  light  of  the 
values  of  rmo„  given  in  Table  1,  it  is  conceivable  that  such  a  modifica¬ 
tion,  even  if  only  a  few  centimeters  in  height,  could  be  extremely 
effective. 


VI.  RECOMMENDATIONS 


One  of  the  limitations  of  the  foregoing  analysis  is  that  physical  and 
thermodynamic  properties  of  the  exhaust  plume  are  taken  to  be  weighted 
averages  of  the  properties  of'  the  constituent  molecular  species. 
Realistically,  the  number  density,  mean  free  path,  and  collision  rates  of 
each  molecular  species  will  not  necessarily  be  equal  to  one  another  at 
any  specified  point  in  the  flow  field. 

Further,  in  this  analysis,  the  effect  on  the  exhaust  plume  of  the 
atmosphere  through  which  the  spacecraft  is  traveling  has  been 
neglected.  Although  the  molecular  number  densities  associated  with  low 
earth  orbit  altitudes  are  small,  the  velocity  of  these  molecules  relative 
to  a  spacecraft  advancing  in  orbit  is  significant.  The  relative  approach 
of  these  atmospheric  molecules  toward  the  spacecraft  is  parallel  to  the 
plane  of  the  nozzle  exit  and  therefore  nearly  perpendicular  to  the  locus 
of  points  describing  the  breakdown  surface  of  the  exhaust  plume. 
Collisions  between  free  stream  molecules  and  exhaust  gas  molecules 
occurring  in  the  transition  region  of  the  exhaust  plume  will  result  in  a 
transfer  of  energy  to  molecules  in  this  region. 

The  resolution  of  each  of  the  limitations  above  will  further  improve 
the  accuracy  of  the  results  presented  in  this  analysis. 


APPENDIX  A 

ON  THE  GRADIENT  IN  CENTERED  RAREFACTION  FAN 


Joseph  Falcovitz 
Feb.  19,  1985 

Purpose:  Evaluate  the  Bird/Knudsen  number 


along  streamlines  of  centered  rarefaction  fan. 


Reference:  J.A.  Owczarek,  "Fundamentals  of  Gas  Dynamics” 

(1)  Mach  wave:  sin  y  =  1/M 

(2)  p  •=  p  ♦  0;  0  =  streamline  angle 

For  a  centered  fan,  flow  field  is  function  of  0. 

Denote:  r  =  distance  from  fan  vertex 
s  =  coordinate  along  streamline 
p  =  density 

pf  =  stagnation  density 

Then  (dp/ds)  is  a  (-sin  y)  component  of  : 

r  d/8 


(3)  Sfi-. 
ds 


1  dp  . 

T^s,n  A*  = 
rd/8 


_1_  dp 
rM  d/? 


Denote  the  non-dimensional  gradient: 


(4,  «(M)  =  S  § 


-•  >  -4  -•  -T  -4 
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Then  for  a  centered  fan: 


(5) 


*(M)=  _  J .*£  =_i. 

Mp  dp  M p  dM 


d£ 

dM 


-i 


da  _ 
dM 

4ft  « 

dM 


A  _  1 

— . —  »  *  ,  2v'/2 

M2cos  M  (1-l/M  ) 

2  J/* 

M[l  +  (7-1)M!/2] 


_ (y+i)M 

3M  =  dM*dw”  2(MJ-1)'/J[1+(7-1)mV2] 


P,  =  p[1+(7-1)m72]'/<,") 


=-pM[l+(7-1)M2/2] 


=  -M[.+(7-1)MV2] 


(6)  k(  M)  =  2Vm2-1/M(7+1) 


APPENDIX  B 

EVALUATION  OF  INTEGRALS 


Given 


(djp/dfi)  =  f  nTcm7r"3/2c3e“(c2+u2_2cu  c©* 

J  o 

Via  the  method  of  Noeller  [Ref.  16],  the  substitution  k  =  c  -  u  cos 
is  made.  The  right-hand  side  of  the  equation  becomes 

_  oo 

nrcmn~V!  f  (k+u  cos  *>  dk 

-  U  COS  iff 


which  in  expanded  form  equals 

,  rtoo 


nTcm7r~3/2e  u2euZ  cos2  k3e“k2  dk  +  3  u  cos  if  f  k2e 

W  _ _ mJ.  —  >  •  n 


-  U  COS  if 


-  U  COS  if 


+  3  u2  cos2  ^  J  ke"k2  dk  +  u5  cos3  i>  J  e~k2  dk 

-  u  cos  i  -  u  cos  i 

The  evaluation  of  these  four  integrals  is  as  follows: 


f°° 

1  •  k3e~k2  dk  =  (u2  cos2  i  e”u2  cosZ  *  +  e“u2  cos2  1')/2 

-  U  COS  if 

pOO 

2.  k2e“k2  dk  =  -(u  cos  i  e“u2  cos2  ^)/2  +  (Vrv'4)erfc(-  u  cos  i) 

-  U  COS  if 

C°° 

3.  I  ke“k2  dk  =  (e"u2  cos2  V)/2 


s 
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LIST  OF  SYMBOLS  IN  PROGRAM  "FLUX" 


Input  Parameters 


Symbol 

Name 

Units 

Description 

GAMMA 

1 

average  specific  heat  ratio 

MW 

M 

kg/kmole 

average  molecular  weight 

DELTA 

6 

m 

average  molecular  diameter 

R1 

m 

spacecraft  radius 

XI 

m 

nozzle  half-width 

ME 

M. 

exit  plane  Mach  number 

DM 

dM 

incremental  Mach  number  step 

MMAX 

Maw 

estimate  of  maximum  flow  field 
Mach  number 

TSTG 

Tt 

K 

stagnation  temperature  of 
exhaust  plume 

DSTG 

P, 

kg/m* 

stagnation  density  of 
exhaust  plume 

P 

P 

breakdown  parameter 

IPHI , 
PHIDIV 

number  of  ^-divisions  for 
integration 

IPTS 

number  of  target  points 

XMIN 

m 

x-coordinate  of  target  point 
nearest  nozzle  lip 

XMAX 

m 

x-coordinate  of  target  point 
furthest  from  nozzle  lip 

Other 

Parameters 

Symbol 

Name 

Units 

Description 

PI 

ir 

3.1415926 

E 

e 

2.7182818 

RU 

universal 
gas  constant 

J/mole  •  K 

8314.41 

NO 

Avogadro's 

number 

molecules 
per  kmole 

6.022045  x  102‘ 
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SIGMA 

a 

m* 

collision  cross  section 

RG 

R 

J/kg.K 

R  =  31/M. 

XDIV, 

LRATlO. 

XCOORti, 

MM 

variables  used 
in  determination 
of  intermediate 
target  points 

Cl 

2/(1  -  1) 

C2 

(I  ♦  1)/(T  -  1) 

C3 

(r  -  i )/(*  ♦  i) 

C4 

(r  -  n/2 

C5 

Or/*)1'* 

C6 

i/(r  -  D 

NU1 

vi 

radians 

Prandtl -Meyer  function  at 
nozzle  exit 

RMAX 

rffl« 

m 

value  of  r  corresponding  to  Mmo. 

PHIMIN 

^mln 

radians 

minimum  value  of  0  used  in 
integration  procedure 

DPHI 

d* 

radians 

incremental  0-step 

PHI 

* 

radians 

spatial  coordinate 

FLUX 

ip 

molecules 
per  m2  per  sec 

flux  arriving  in  vicinity  of 
target  point 

ML 

M 

local  flow  Mach  number 

NU 

V 

radians 

Prandtl-Meyer  function 

THETA 

8 

radians 

angle  between  spacecraft 
surface  and  flow  streamline 

BETA 

& 

radians 

angle  between  spacecraft 
surface  and  Prandtl-Meyer  wave 

R 

r 

m 

distance  from  nozzle  lip 

DR 

dr 

m 

incremental  distance  step 

MPS 

cm 

m/sec 

most  probable  thermal  molecular 
speed;  (2RT)1^ 

UMAG 

\u\ 

magnitude  of  stream 
flow  velocity 

LDOTU 

LMAG 

LNDS 

.L*  u 

L-  ndo 

variables  defined  by 
equations  3.33,  3.34, 

3/35,  and  3.43. 

PSI 

0 

radians 

flux  back  flow  angle 

DOMEGA 

dfl 

steradians 

solid  angle  subtended  by 
portion  of  breakdown  surface 

u 

u 

^unenjj^inless  ^ow  speed. 

C7 

u  cos  f 

C8 

u*  cos2  # 

C9 

CIO 

eU2  cos2  # 

NO 

"T 

molecules/m1 

number  density 

DFLUX 

dip 

molecules 
per  m2  per  sec 

incremental  flux  directed 
toward  target  point 

TFLUX 

JP 

molecules  total  flux  incident  near 

per  cm2  per  sec  target  point 

MLIM 

Mimi 

maximum  stream  Mach  number 
for  which  breakdown 
is  significant 

uu 
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